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Abstract 

The Wilson Coefficients for all 4-parton operators which arise in matching QCD to 
Soft-Collinear Effective Theory (SCET) are computed at 1-loop. Any dijet observable 
calculated in SCET beyond leading order will require these results. The Wilson coef- 
ficients are separated by spin and color, although most applications will involve only 
the spin- averaged hard functions. The anomalous dimensions for the Wilson coefficients 
are given to 2-loop order, and the renormalization group equations are solved explicitly. 
This will allow for analytical resummation of dijet observables to next-to-next-to- leading 
logarithmic accuracy. For each channel, there is a natural basis in which the evolution 
is diagonal in color space. The same basis also diagonalizes the color evolution for the 
soft function. Even though soft functions required for SCET calculations are observable 
dependent, it is shown that their renormalization group evolution is almost completely 
determined by a universal structure. With these results, it will be possible to calculate 
hadronic event shapes or other dijet observables to next-to-leading order with next-to- 
next-to-leading log resummation. 



1 Introduction 



While the main goal of the Large Hadron Collider at CERN is to find evidence for physics 
beyond the standard model, the standard model itself must be accurately understood before 
any new physics claim can be made. The most common way of comparing data to predictions 
of the standard model is through Monte Carlo simulations. While these simulations are 
remarkably accurate in reproducing the gross kinematic features (and often the fine structure) 
of observables, they are formally accurate to only leading-order in perturbation theory and 
with resummation of only the leading Sudakov double- logarithmic singularities (LL). While 
much theoretical work has been devoted to computing phenomenologically relevant observables 
to next-to-leading order (and sometimes even next-to-next-to-leading order), very little is 
known about the effects of higher-order resummation. Since observables at the LHC generally 
involve many scales, it is likely that resummation will lead to a quantitatively significant 
improvement of our theoretical predictions. Recent advances in Soft-Collinear Effective Theory 
have shown resummation to be achievable at the next-to-leading (NLL) [IJE], next-to-next- 
to-leading (NNLL) [31 IU El E] and next-to-next-to-next-to-leading (N 3 LL) logarithmic order 
[3 El El QUI ttU C21 [13] . Practical results have included a precise measurement of a s from event 
shapes [TOl ITT], better agreement with Tevatron data and reduced theoretical uncertainty for 
the the direct photon pt spectrum [5], and a bound on possible new colored states [TJ]. In 
all of these cases, higher order resummation was critical to the improved theoretical precision, 
and so it is natural to ask whether resummation might be similarly important for pure QCD 
events. 

In this paper, we take the first step towards the calculation of dijet observables at next- 
to-leading order (NLO) with NNLL resummation. The computation of a dijet observable in 
SCET requires a factorization formula which convolutes parton distribution functions (PDFs), 
Wilson coefficients, jet functions and a soft function. Although the jet and soft functions are 
observable dependent, the Wilson coefficients are universal and independent of the details of 
the factorization theorem. For an NLO calculation, the complete set of NLO Wilson coef- 
ficients for all the 2—7-2 processes of QCD is required, but explicit expressions have only 
been presented in the literature for the qq — > qq channel and its various crossings [15]. The 
necessary diagrams for qq — > gg and its crossings were computed in [T6l ITT] and for gg — > gg 
in [TTJ and [T8]. These amplitudes can be used to extract the Wilson coefficients for the 
remaining channels. This computation, and the compilation and simplification of all of the 
1-loop expressions, is a main result of this paper. 

The other main result of this paper is the elucidation of how the color structures and spin 
states interfere in the context of resummation. Processes with four hard partons in QCD 
separate into a number of color structures, corresponding to separate operators in SCET, 
which mix under RG evolution. Similar color mixing has already been studied for collinearly- 
regulated soft Wilson lines using the traditional approach to resummation [19], and for quark 
scattering [15], gauge boson production [20] and tt production [6] using SCET. For NNLL 
resummation, the renormalization group equations must be solved to least 2-loop order. A 
general form for these RGEs is known [201 EH E21 ESI EH EH] , as is the related soft function 
evolution equation |19j . 
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Since jet functions and PDFs are color diagonal, the color mixing in the evolution of the 
Wilson coefficients is exactly compensated for by color mixing in the soft function evolution. In 
fact, since the Wilson coefficients are universal (observable independent), the mixing terms in 
the soft function must also be universal. Moreover, all of the color-mixing is proportional to the 
universal cusp anomalous dimension (up to at least 2-loop order). Therefore, the color-mixing 
effects can be diagonalized once-and-for-all by a careful choice of operator basis. This natural 
basis for evolution is different from the one that is natural for matching, in which the Wilson 
coefficients take a particularly simple form. This implies that although the color mixing in the 
evolution can be diagonalized, there are still additional interference effects relevant beginning 
at NNLL due to the non-diagonal components of the NLO hard and soft functions. 

This paper is organized as follows. In Section [2} we give an overview of factorization for 
dijet processes in SCET. Section [3] introduces the operators for the 2 — > 2 processes. Their 
Wilson coefficients are computed in Section H] and the RGEs are presented and solved in 
Section [51 Then, in Section [6] we explain which features of the soft function evolution are 
universal. Finally, [7| presents some conclusions and possible applications of these results. 

2 Factorization in SCET 

Cross sections for dijet observables calculated with Soft-Collinear Effective Theory [29 | 1301 13T] 
will have the general form 

da-dn -^—cy-s JI -c l I ®j®j®f®f. (1) 

i,J,r, mit 
channels 

Here, dll is the Lorentz-invariant phase space and N init are the number of initial states, which 
are averaged over in computing the cross section. Cj are Wilson coefficients, with / indexing 
the different color structures and V indexing the different spins. The Wilson coefficients encode 
information about the hard 2 — ^2 scattering process. Sjj are the soft functions. These form 
a matrix in color space, but are independent of spin. The jet functions J and the PDFs / 
are spin-averaged and color- independent. The sum is also over channels, which for 2 — > 2 
processes are qq — > qq, qq — > gg and gg — > gg and their various crossings (like qg — > qg). 
Everything has an implicit channel index, which we suppress for clarity. One can show that 
at leading power in the SCET expansion the channels do not interfere and the cross section 
for each channel can simply be added together. For more details, see for example [5]. 

The effective field theory calculation begins by matching from QCD to SCET. This involves 
enumerating the appropriate operators in both theories and calculating matrix elements in 
relevant final states. Wilson coefficients for the SCET operators are then adjusted so that 
they reproduce the amplitudes from QCD to the desired order in perturbation theory, if 
evaluated at the hard scale /i/j. If the SCET operators are evaluated at lower scales, their 
matrix elements change, and include the effects of resummation. Matrix elements of outgoing 
collinear fields give jet functions, incoming collinear fields give PDFs, and matrix elements 
of soft Wilson lines give the soft function. Due to factorization, all of these objects can be 
evaluated independently and combined together using a factorization formula. 
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To illustrate the procedure in detail, consider the process qq' — >■ qq' . We will label the 
incoming partons as 1 and 2 and the outgoing partons as 3 and 4. The operators in QCD 
appropriate for this channel have the form 

0?r CD = {mi^q 2 ){mi^' qi ) . (2) 

Here, / indexes the color structure (T\ = r a or T 2 = 1), and T and V index the spin (e.g. 
r = T' = Pjr, = ^ (1 — 75))- The operators in SCET follow from replacing each QCD quark 
with a collinear quark and adding soft and collinear Wilson lines, 

Cfr ET = (x^Ta.TY^i^Ta^xi) • (3) 

The fields denoted by x% are collinear quarks in the raf direction wrapped in rij-collinear Wilson 
lines (not shown). The are soft Wilson lines in the fundamental representation of SU(3). 
For more details of the notation, see for example [5]. 

In SCET there is a separate set of collinear degrees of freedom associated with each di- 
rection 1234, and a single set of soft degrees of freedom which can communicate among the 
collinear sectors. When we calculate a cross section from the squared matrix elements of 
SCET operators, contributions from different collinear sectors do not interfere and can be 
factorized into separate calculations. Vacuum matrix elements of collinear fields associated 
with outgoing directions turn into jet functions. For example, we may be interested in the 
inclusive jet function in the 4 direction. This requires the evaluation of 



<0|x4(x)r X4 (0)|0>=tr 
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Gluon jet functions are defined similarly [5], and less inclusive jet functions are possible as 
well [321 E31 131] • When the collinear field is associated with an incoming direction, the matrix 
elements should be taken in the state of the appropriate nucleon (Ni or N 2 ), which produces 
parton distribution functions (PDFs). For example, the incoming collinear field in the 1 
direction gives 

(^|Xi(a?+)rxi(0)|JV 1 > = ifii-Pitr[ji 1 ri /^/^(Oe 44 ^^. (5) 
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Here, fq/^iO is a PDF. The dependence on x + = ri\ -x^- indicates that, in contrast to the 
jet functions, fluctuations of initial states in the transverse direction are power suppressed. 

The final contribution to a cross section from the SCET operators comes from the soft 
Wilson lines. For the operator in Eq. fl3]), these are 

W, = T I (y£TiY 2 )\ (Y 3 t T/Fi)* 3 ii I 

= T{{YlT I Y 2 ){Y^T I Y l )) . (6) 
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Wj has 4 fundamental color indices which have been suppressed in the second line since the 
index structure is clear. Taking vacuum matrix elements of these Wilson lines produces a soft 
function 

Su ({*}, <) = J2 (°l ^ \ x ») (*s\ >Vj |0) F s ({k}) , (7) 

x 3 

where the sum is over soft radiation in the final state. The function Fs({k}) encodes the 
dependence on various projections related to the definition of the observable. At tree level, 
the soft wilson lines are trivial (1$ = 1) and the function Fs({k}) must reduce to a product of 
delta functions since there is no soft radiation. For the example being considered (qq' —> qq'), 
we have Wi = r a and W2 = 1 and the tree level soft function takes the form 



2 u 

c\ 



l[S{k). 

{k} 



Note that the jet functions and PDFs are spin averaged, and the soft function is indepen- 
dent of spin. Thus, the spin correlations are determined completely at the hard scale and are 
stored in the Wilson coefficients. To be explicit, the SCET Lagrangian contains a sum over 
all the possible spin states T, 

C = J2 C i°i- ( 9 ) 

When one calculates a matrix-element squared from this Lagrangian, the operator products 
turn into jet and soft functions which have no spin dependence. So, we can write heuristically 
that 

\ M I 2 = w~ E c ^ ■Sji-c ] i®J®J®f®f- (10) 

iVinit J,j,r 

This is the same as Eq. (CE]) once the phase space factor is added. Therefore, for a spin-summed 
cross section we only need the spin-summed hard function defined by 

H IJ = Y J C V I C T /. (11) 
r 

Then, 

\M\ 2 = ^-Y< HijSji ® J ® J ®f®f- ( 12 ) 

7Vinit 1, j 

Since the spin information is retained in the Wilson coefficients, spin correlations can be 
studied using SCET, if desired, simply by summing only over the desired spins T. 

As a check on the Wilson coefficients, we can compare the 2—7-2 scattering cross section in 
SCET to the spin-summed cross section in full QCD. If the matching has been done correctly, 
these two calculations should agree. Although the 2—7-2 processes are infrared divergent 
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and not physical, the comparison can be made in 4 — 2e dimensions as a formal check. Since 
there is no structure to the outgoing or incoming partons for 2 — > 2 parton-level scattering, 
only virtual graphs contribute. Thus, the jet functions, soft functions, and PDFs are fixed to 
their tree-level values, which produce delta functions of the parton momenta. Only the soft 
function provides non-trivial structure, through its color factors. Putting these observations 
together, we should find 



m qcd 



init 



Err otrcc 
H U b JI 



I.J 



Ni 



-Tr[HS 



trccn 



(13) 



init 



We will use this equation as an important check on our calculations, comparing the matched 
SCET prediction to the cross sections from Ellis and Sexton [28J. 



3 Operators 

In this section, we present the operators relevant for dijet production. Again, we will label 
the incoming partons as 1 and 2 and the outgoing partons as 3 and 4. 

First consider the channels involving only quarks and anti-quarks. We start with qq' — > qq', 
with q and q' different quark flavors, e.g. u(pi) + d(j>2) — > u(ps) + d(p$). In this channel, there 
are two separate color singlet operators for each spin 

O^v = (x^y^TY^ixsY^r^r'Y.Xi) (14) 
= (X4Yh^Y 2 X2)(X3Yh^'Y lXl ) . (15) 

Here, r a refer to generators of SU(3) and the color adjoint index a is summed over. V and V 
label the spin of the operator. 

Since QCD is non-chiral, left- and right-handed fermions can be thought of as separate 
species; thus it is simplest to take a basis where Y and V are either P L = | (1 — 75) or 
Pr — \ (1 + 75); as in [15]. Note that the primes on the T go with the second contraction, not 
with the second quark species. For example, for the u(pi) + d{p2) — > u(pz) + d{jp^) process, 

{LL = P L P L ] ( u L d L -> u L d L ] 

LR = P L P R \ refers to I u R d L -> u R d L \ . (16) 

and so on. 

For the crossed processes the fermions are contracted differently. For example, in qq' — > q'q, 
the fields labels change from 4231 to 3241. The new operators are 

0{f T , = txzYy^TY^iXAYlT^T'Y^i) (17) 
Olrr' = {xzYh^Y 2X2 ){xAYlrV'Y lX i) . (18) 
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There are 4! possible permutations of the labels 1234, but the Wilson coefficients for all chan- 
nels are determined by crossing relations of the Mandelstam invariants s, t and u. Therefore, 
we label the operators with the appropriate permutation of stu. The momentum routing and 
operators for various crossings are given in Table [TJ Our conventions are chosen to agree 
with [28]. Note that the stu channel is defined as qq' —> qq' which has a t-channel pole at 
leading order. 

As long as the quarks are distinguishable, only one set of operators is required for each 
channel. However, there are additional subtleties when identical particles are present. Con- 
sider the process qq — > qq. At tree-level, this process gets a contribution from the t-channel 
and the li-channel. The t-channel is like qq' — > qq' and can be matched onto Of^, while the 
M-channel is like qq' — > q'q, which matches onto Offi. As long as there is no interference, both 
of these operators can be turned on at the same time in the Lagrangian, and the amplitudes 
added incoherently. In fact, for the LR and RL operators, in which the quark helicities are 
different, this is the case. However, for the LL and RR, interference is important. This can 
be seen in the effective theory because the operators are not linearly independent. As was 
shown in [T5j . 

OTI^BuOTll, (19) 

where 

B LJ = I 




This crossing matrix satisfies B 2 = 1, and includes the —1 from Fermi statistics. Thus, to 
treat the LL and RR channels, only the Of^ L operators are required. We will use this matrix 
to extract a single set of Wilson coefficients for the qq — > qq channel in the next section. 

For the gg — > qq processes, there are three color structures and 16 independent spin/helicity 
combinations. For the quark spins, it makes sense to project out the left- and right-handed 
states using the r = Pl or Pr projectors. For the gluon helicities, we could use a similar 
projector formalism, but is is easiest to simply write ± as a label for the (incoming) gluon 
helicities. For the gg — > qq channel, the operators in QCD are then 

Og ± = (q 3 T a T b Tq 4 )(A a l ± A b ±) 

O s £ ± = (g 3 rV a rg 4 )(^ ± ^ ± ) (21) 
O s 3 t ± = (^ ab Tq 4 )(A a l ± A b 2 ± ), 

<D% ± = (q 3 T ab rq i )(AfA b± ) , (22) 



which we generically write as 



with Ti = r a T b , T 2 = T b T a and T 3 = 5 ab . The T dependence of these operators determines the 
spin of the out going quarks, and since QCD is non-chiral, the amplitude will be independent 
of T and so we leave the T index on the operators implicit. The following equation should 
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make the helicity conventions clear. 



±± 



L = P L ; 
L = P L ; 
R = Pr] 



> refers to < 



g+g^ -> q L q R 
g + g~ q L q R 
g + g- q R q L 



(23) 



It is straightforward to translate these QCD operators into operators in SCET. Gluon 
fields with associated collinear Wilson lines are labeled A^ ± or equivalently A^ ± . To make the 
notation clear, A^ annihilates an incoming positive (+) helicity gluon and creates an outgoing 
negative (-) helicity one. We write y^ b for a soft Wilson line in the adjoint representation in 
the direction n M , as in [30]. Then, performing a soft-field redefinition 



leads to 



U I±± 



xlA^TA^xi) (Y.]yr'Tf b 'y^ b 'Y i 



(24) 



(25) 



The indices a and b are adjoint color indices and i and j are fundamental color indices. We 
have also used factorization to pull the soft Wilson lines away from the collinear fields. When 
these operators are squared, the a, b, i and j indices in the collinear part (left brackets) are 
contracted with their counterparts in the adjoint operator through delta functions, since the 
collinear interactions are color diagonal. Therefore, the a, b, i and j indices in the soft part 
(right brackets) also get contracted with the soft part of the adjoint operator. Thus, collinear 
fields can be completely ignored when calculating color-related effects. This situation is almost 
identical to the direct photon case, except now with multiple color structures; more details 
of how such a factorization arises and how the color and spin indices are contracted can be 
found in [5]. 

The crossed operators are constructed in a straightforward way by permuting the indices. 
For example, for qg — >• qg, the operators are 



(26) 



Various crossings and the momentum routing conventions, which again are chosen to agree 
with (28] , are shown in Table |2j 

For gg — > gg, we also label the operators by the gluon helicities, rather than putting in 
explicit helicity projectors. There are 8 color channels for this process, 16 spin states, and no 
crossings. The matching is more conveniently performed using an over-complete basis of the 
following 9 operators 



/H±±;±± / /1±Q /l±6 /l±c A±d,\ I rp a' 'b' c' 'dl -\-taa' -\ibb' -\-icc' -\-idd' 

Uj = [A 1± a 2± a 3± a 4± j yi! y l y 2 y 3 y A 



(27) 



where the color structures are given by 

rpabcd _ rp^^fycydj 



rpabcd Ti_r_<i_c_6_dl 



TrfrVVV 
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rpabcd 
1 2 


= Tr| 


"^-i ^-b r-d ^cl 
T T T T \ 


rpabcd 
1 7 


= Tr 


[r<v; 


]Tr [r c T b \ 


rpabcd 
1 3 


= Tr| 


T T T T \ 


rpabcd 
1 8 


= Tr 


[tV; 


Tr [r c T d 


rpabcd 
1 4 


= Tr| 


T T T T \ 


rpabcd 
1 9 


= Tr 




Tr[rV 6 ; 


rpabcd 
J 5 


= Tr| 


'^a—C^d^bl 
T T T T \ . 











While this basis can be reduced to a linearly independent set of 8 operators, there is no need 
to do so. As long as the interference effects among the 9 operators are properly accounted for, 
there is no difficulty in the effective theory with the basis being overcomplete. 

4 Wilson Coefficients 

The Wilson coefficients Ci for the operators Oi are calculated in matching from QCD to 
SCET. They are uniquely fixed by requiring that the SCET Lagrangian reproduce the matrix 
elements of QCD order-by-order in a s . Since virtual graphs in SCET are scaleless, they vanish 
in dimensional regularization. Moreover, since infrared divergences in SCET are the same as in 
QCD, they cancel in the matching step. Therefore, the matching coefficients can be extracted 
entirely from IR-regulated virtual graphs in QCD. All of the following results are derived using 
dimensional regularization and the MS subtraction scheme. 

There are a number of results in the literature which are relevant for the calculation of 
the Wilson coefficients. The classic work [28] presents the cross sections for all the channels 
separately, summed over spins and colors. This is not enough information to extract the 
Wilson coefficients, which we need at the amplitude level. Refs. [36] and [17] have calculated 
the necessary diagrams for the qq — > qq channels. Refs. [IB] and JT7J have the qq — > gg results. 
For gg — > gg most of the loop amplitudes appear in [T7J, except for a few of the color channels 
which can be found in [18]. Much of the previous work is written using the spinor helicity 
formalism. In order to simplify the Wilson coefficients, we will convert the spinor helicity 
products into Mandelstam invariants. 

Our momentum convention is p\ + p 2 — > pz + Pi- We define Mandelstam invariants as 

s = (Pi + P2)\ t = (pi - p 3 ) 2 , u = (p 2 - p 3 ) 2 . (29) 

This convention agrees with [28], but differs from Chiu et al. [15] via u -H- 1. For all processes, 
physical kinematics has s > and t,u < 0. Since the matching is done at the hard scale we 
will indiscriminately use 

s + t + u = (30) 

to simplify the Wilson coefficients. 
4.1 qq — > qq channels 

First, consider the 4-quark processes. There are six channels involving two different flavor 
quarks (such as u and d, which we call generically q and q') 

qq' — > qq', qq' —> qq', qq — > q'q', qq' — > q'q, qq — > q'q, qq — > q'q' . (31) 
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12 34 


crossing 


operators 


12 ->> 34 


crossing 


operators 


yy r yy 


stu 


vy4^ 1 y2Ay3 ± ± Hi) 


yy ' y y 


sut 




qq' — > gg' 


uts 


(g 2 TTg 4 )(g 3 r / Tgi) 


qq' —> q'q 


tus 


(q 2 TTq 3 )(q 4 T'Tq 1 ) 


gg — > g'g' 


tsu 


(g4Trg 3 )(g 2 r / Tg 1 ) 


qq — > q'q' 


ust 


(q 3 TTq i )(q 2 T'Tq 1 ) 



Table 1: Crossing relations for the 4-quark channels. For example, the Wilson coefficients for 
the tsu channel are obtained from the stu channel via Cf. tsu (s,t,u) = Cf. stu (t, s,u). 



The first channel we call stu. The others are related by crossing symmetries, as shown in 
Table [TJ In addition, there are three identical particle channels 



qq — > qq, qq — > qq, qq — >■ qq . 



(32) 



All other channels, such as qq' — > qq' and qq — > qq, are the same as one of these by charge 
conjugation invariance of QCD. Charge conjugation also swaps L <H- R, although this is 
unimportant since L and R states never interfere. 

For all the non-identical particle channels, the Wilson coefficients can be written as in the 
general form 





\s,t,u) 




1 + 


4:71 


-2C F L(t) 2 + X x {s, t, u)L(t) + Y 4 


v 2 


C L X R 


[s,t,u) 


C,U 

= 2 " 7 ' 


!- 


a s 
Air 


-2C F L{tf + Xi(s, t, u)L(t) + Y H 


-(2C F 


C L 2 L 


[s,t,u) 






X 2 ( 


C 

s,t,u)L(t) 2C Z(s,t,u) 


}■ 






[s,t,u) 




k 4vr 


x 2 


C 

s,t,u)L(t) + 2C ^Z(u,t, s) 


}■ 





with 



Xx(s,t,u) = QC F -Pq + 8C f [L(s) - L(u)] - 2C A [2L(s) - L(t) - L(u)] 
2C 

X 2 (s,t,u) = —^[L(s)-L(u)} 



C A 



Y = C a [ 1 { + A+cJ^-16 



Z(s, t,u) = -( [L(u) - L{t)f + 2[L(u) - Lit)} + ^1+^ 

s \ s s 



(33) 
(34) 
(35) 
(36) 



where 



Po = y Ca ~ 3^/ 



(37) 



and Ca = N = 3 and Cp = N 2N l = §, with N = 3 the number of colors. The function L(x) is 
defined by L(x) = log(|x|/;U 2 ) — m9(x). For massless 2—^2 scattering, s > and t, u < 0, so 



L(t) 


= log 


-t 

/I 2 




L(u) 


= log 


— M 
/i 2 


(3 


L(s) 


= log 


s 

— -nr. 
fi 2 





It is important to cross the Wilson coefficients before extracting the imaginary parts from 
these logarithms. 

The hard function Hpj = Ylr^jCj* * s then, up to order g 4 a s , 



H n (s,t,u) = 8g 



4 s z + u< 



t 2 



2 4 a s 



S Z + u- 



(-4C F L(t) 2 + 2X x (s, t, u)L(t) + 2Y \ (39) 



s 2 ?/ 2 
+ -(CU - 4C F )Z( S ,t,n) - -(2CU - 4C F )Z(n,t, S ) 



^ 13 ( a ,t ) u) = 8^ 



+ 



4tt L t 2 
H 21 (s,t,u) = H 12 (s,t,u) 
H 22 (s,t,u) = 0. 



X 2 (s,t,u)L(t) 



s 2 C 



Z(s,t,u) + ——-Z(u,t,s) 



t 2 2C A 



t 2 2C A 



As a check, we can compare |-M(2^2)| 2 = jf— Tr[HS tree ] to the 2—7-2 matrix-elements- 



squared in QCD, as in Eq. ([TBI . For example, in the qq — > q'q' channel (tsu), at tree-level 

1 ' -LR _ nRL ~ oU 



C{ L =C? R = 2g 2 s -, C^ K = C? L = 2g 



C 2 r = 



(qq -> q'q' 



(40) 



So, 



rrtrcc 



10Wi 2 + !l 2 





At tree-level, the soft function is 



(qq -»■ gV 



(41) 



rrtrec 
7J 



\CaCp 
c 2 



(42) 



For this channel, A^ = 4A^ 2 giving a final result 



M 



(qq^-q'q') \ 



4AT2 



t 2 + u 2 



+ ■ ■ 



(43) 
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The result is the correct tree-level matrix-element-squared for parton level qq — > q'q' scattering 
in QCD. Using the NLO hard function, the result agrees with [25] . 

For identical particles, the Wilson coefficients should be included only for linearly inde- 
pendent operators. This is not strictly necessary, but avoids having to compute interference 
effects. As discussed in Section [31 for qq — > qq in the LL spin configuration, the operators 
Of£ L and Of^ L are not independent. To correct for this, we can match onto the operators 
separately, and then use Eq. (JTHJ) to write 

C LLLL = C RRRR = Q LL ^ ^ u) + B ^L^ % f) {gg gg) {U) 

We have added more spin labels to make the counting of independent states more transparent. 
We can also cross the RL channel into the stu basis so that it can be evolved with the uts 
RG kernel (see Section [5TT1 below) . So we define for the qq — > qq channel 

C LRLR = qRLRL = C LR^ ^ ^ ( gg ^ gq) ? (45) 

c lrrl = qRLLR = Bl jC L j R (s, u, t) {qq -> qq) . (46) 

These contributions can then be combined incoherently. Summing over the six contributions, 
the tree-level hard function is then 

Jg^_ (* + C? A u* + a\t-C A uf -C F {t* + s\t* - C A tu))\ 
IJ ~ C\W\ -C F (t 4 + s 2 (t 2 -C A tu)) C 2 F t 2 (s 2 + t 2 ) )' m m) 

Multiplying by the tree-level soft function and dividing by AN 2 to average over initial states 
gives the correct tree- level qq — > qq matrix element in QCD 

4 ( s 2 + u 2 t 2 + u 2 \ 8 u 2 



\M^„r=- g {—+—)-^ Vt +--- m 

Using the NLO hard function, \Ai\ 2 also agrees with the NLO result in full QCD [28] . 

For qq — > qq, with the momentum convention in Table [U we need the uts and ust operators. 
As for qq — > qq the LR spins add incoherently and should be treated as separate channels, 
while the LL operators interfere. In this case, the same matrix applies. Thus, we need 

(48) 
(49) 
(50) 

These four contributions can then be combined incoherently into a single hard function, as for 
qq — > qq. For qq — > qq, the channels are tsu and tus. 
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Table 2: Crossing relations for the gg — > qq channels. 



4.2 — > qq channels 

Next, we give the Wilson coefficients for the gg — > qq channel and its crossings. There are 
again six crossings 

gg^qq, qg-^qg, qg-^qg, gg -> qq, qg ^ gq, qg-^gq- (51) 

The relevant operators, momentum conventions, and crossing relations are given in Table [2J 
The Wilson coefficients for qq — > gg are the same as gg — > qq, but one must keep in mind 
the different number of color and spins in the initial state for a spin-averaged cross section. 
The amplitudes for quarks are the same as the amplitudes for anti-quarks, but with opposite 
gluon helicities, + -H- — and opposite spins L •H- R. Since we will be mostly interested in 
the hard function Hjj which sums over helicities, the convention we choose for ± and L/R is 
unimportant. 

The Wilson coefficients Cj ' = C XlX2,X3Xi have 4 helicity labels A„ and one color label, /. 
Also, helicity is conserved for massless quarks along a fermion line, so the Wilson coefficients 
vanish unless A4 = —A3. Parity invariance implies that 

£7AiA2;A3A4 q— Ai— A2;— A3— A4 (52) 

Therefore, we only give expressions for half of the Wilson coefficients C XlXo ~ = C XlXr ' + ~ for 
gg ~^ qq with fixed quark helicities A 3 = + and A 4 = — . The other half is obtained through 
Eq. <|S2J». 

The Wilson coefficients for the 4 independent spins and 3 independent colors are 



Cr + (s, *, u) = 2g 2 ^ {l + ^ [-{C A + C F )L{sf + V 1 {s 1 1, u)L(s) + W l {s 1 1, u)] } (53) 
C+- (s, t,u) = 2g 2 - s ^{l + ^- [-(C A + C F )L{s) 2 + V x (s, t, u)L(s) + W 2 (s, t, 
C+ + (s,t,u) = Cr(s,t,u) = 2g\^^W 3 (s,t,u) 
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With these Wilson coefficients, it is easy to compute that the hard function for gg — > qq is 
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and E n = E u . 

The tree-level soft function in this color basis is 
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Using these results, we can check the 2—7-2 cross section. For gg — > qq, 
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which is the correct tree level result. The 1-loop hard function when combined with tree- level 
soft function reproduces the NLO cross section in [28] . 



4.3 gg — > gg channel 

Finally, we give the Wilson coefficients for gg — > gg. As discussed in Section [31 there are 8 in- 
dependent color structures, but we match to an overcomplete basis of 9 color structures. There 
are 16 possible helicity amplitudes for each color structure, giving 144 matching coefficients. 
Fortunately, there is a great deal of symmetry for this process which relates the various helicity 
and color subamplitudes. Parity allows us to give the results for half of the helicities, with the 
other half obtained via Eq. (|52|) . In this section, we will number the helicities T = 1 • • • 16, 
with the correspondence given in Table [3j 

Only the first 6 color structures in Eq. (f28j) are non-zero at tree-level, and then only for 
6 of the possible 16 helicity amplitudes. The tree-level amplitudes, Ai^ are given in Table HI 
At one loop, all 16 helicity amplitudes are non-zero for all color channels. The NLO matching 
coefficients are given by the following formula 
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Table 3: The 16 helicity structures (T) for the gg — > gg channel. 
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Table 4: Tree-level matching coefficients, A4^, for the gg — > gg channel. 



where the expressions for are given in Table [51 The quantities Q] are written in terms of 
A, B, and J 7 , which are as follows: 

A(s, t, u) = -2C A L(u) 2 + ( - 2C A [L(s) - L(u)} + A,) L(u) + - ^) C A + ™n f 



ft?[L(ti) - L(s)} - {[L(u) - L(s)} 2 + vr : 
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s — u 
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1 /s 2 s 2 2s 2s 
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(61) 



With these Wilson coefficients, we compute the hard function for gg — > gg is 
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Table 5: The expressions for for the indicated helicity structures. 
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The 1-loop hard function, Hj^ lo ^ is complicated and is therefore not written explicitly. 
The tree-level soft function in this color basis is 
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with 



a = 3 + C\ b = 3-C 2 A c = C A -3C 2 A + 3 

d = 2C 2 A C F e = C A . (64) 

Combining Hjj with Sj r / e and using N = 3 gives 

I ka ,2 1 „ .tree & A C F <f (s 2 + t 2 + M 2 ) ( g 4 + t 4 + M 4 ) 

l*W«ol = 4^# H "Sji = —32 + " ' ' (65) 

which is the correct tree-level matrix-element-squared for gg gg. As with the other chan- 
nels, the 1-loop hard function for gg gg, when combined with the tree- level soft function, 
reproduces the NLO cross section in full QCD [28] . 



5 Renormalization Group Evolution 

The renormalization group evolution of Wilson coefficients in SCET is determined by the same 
virtual graphs as in the matching step. However, the evolution is not uniquely fixed given only 
the Wilson coefficients themselves; one needs also the mixing terms. The Wilson coefficients 
satisfy a renormalization group equation of the form [22l |23| [2] 



C T M = ») , (66) 



dln/i 

where, up to at least 2-loop order, Tfj has the form 

rfj(s, t, u, y) = ( 7cus P ^- In — ^ + 1h - ^ " J Su + 7cus P M/j(s, t, u) . (67) 
V z fx a s j 

The /3(a s )/a s term compensates for the g 2 (y) dependence in the leading order Wilson coeffi- 
cients. The QCD beta function is (3(a) = —2a (f-) (3 + ■ ■ • with (3 given in Eq. (137|) . The 
cusp anomalous dimension is 7 cusp = 4 f j^J + • • • . For NNLL resummation, one needs the 
/3-function and 7 cusp to 3-loop order, which can be found in [5], and the 7# and Mjj terms to 
2-loop order, which we give below. 

The specific form of the hard function RGE in Eq. (1671) implies that 



1. All of the fi dependence is proportional to the identity matrix in color space. It is fixed 
by a single Casimir Ch which is known exactly. Explicitly, 



ch = C Ri = n q C F + n g C A , (68) 

i 

where C Ri is the quadratic Casimir of representation and n q (n g ) the total number of 
quarks (gluons) in the hard process. The sum is over both initial and final states, 
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2. All of the color mixing terms are proportional to the universal cusp anomalous dimension 
Tcusp- They are fixed by a single matrix Mu(s,t,u) which depends on the channel and 
is known exactly Once a color basis is chosen, this matrix can be computed from 



M 



T, • T\ 



Lis 



L(t) 



(69) 



where si 2 = s 34 = s, si 3 = s 2 4 = t and s 24 = S13 = it and L(x) is given in Eq. fl38|) . The 
operator Tj • Tj acts in color space on the operator basis, as discussed in [20] . 

7ij depends on the number of quarks and gluons involved, but not on the color structure. 
It is given by 



1h = n q ^ q + n g -f g , 

where, lq = ) (-3C F ) + ■■■ and lg = (%) (-/3 ) + 
anomalous dimensions can be found up to 3- loops in [22] . 



(70) 

. These quark and gluon 



The form of Eq. (ISTI) is known to hold up to 2-loop order. Here, we verify it to 1-loop by direct 
calculation. The 2-loop results are then used to solve the RGE of the Wilson coefficients for 
NNLL resummation. Eq. ( 1671) may hold to more than 2-loops, however at 3-loop order new 
possible terms may arise which appear not to be forbidden by general arguments [21] . 

The universality of the off-diagonal terms in Mjj allows the mixing to be diagonalized once 
and for all. This diagonalization was shown at next-to-leading order in the traditional approach 
in [19], and we find an eigensystem in SCET which is essentially identical. Differences in the 
eigenvalues are due to different conventions for the diagonal terms, related to the way double- 
counting of soft-collinear divergences in the two theories are handled. Also, the explicit factor 
of log ^| in Eq. (!67j) is a convention. We could have used any other scale in this logarithm, 
such as log the physical prediction is independent of our convention. 

Let us denote the linear combination of Wilson coefficients which are eigenvectors of Mjj 
by C l K , with eigenvalues Xk- In the diagonal basis, the evolution equation takes the form 



dln/i 



(A*) 



°JL\ 

7cusp 2 m ^ 



7 + lH + 7cuspA^ 

a. 



(71) 



As we can see, the evolution is now local in color space. Since the off-diagonal mixing terms 
are proportional to 7 CU s P , this holds to at least order a 2 s . Thus, it can be solved in closed form 
once-and-for-all [351 [36] . The solution to the RGE is 



CM = — TT exp 



c H S(/j, h ,/j,) - A H (fi h ,ij,) - A v {n h ,ii) ( Xk + ~~ hi — ^ 

1 



CkM, (72) 



where the functions S*(z/, /i) and A(u,fi) are the same as in [5]. They have the general form 

da' 



S(u,fi) = — [ da ^°" sp ^\ , 

Ja a (u) /3(a) J* s {v) P{<x') 



(73) 
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A r (v, At) 



aM , TcuspH 



(74) 



Ajj(u,fj) is the same as A r (z/, /i) but with 7^ replacing 7 cusp . Closed-form expressions for 
these functions in renormalization-improved perturbation theory can be found in [57]. 

Since the RGE is independent of spin, we can also solve directly for the evolution of the 
spin-summed hard functions Hjj = CjCj*. If we go to the basis in which Mu is diagonal, 
then the hard functions evolution equation is particularly simple 



dln/i 



7cus P c H In 



+ A* + At, + 2 lu - 



2/3(a s ) 



H KK ,(iS), (75) 



where we have used 7^ = jjj. The solution, reinstating the stu dependence for clarity, is 

2c H S (fi h , n) - 2A H (n h , fj) 



H KK ,{s,t,u,ti) = — , 9 exp 



x exp 



a s (jd) 2 

A v {lJ,h,lj) [ X K (s,t,u) + A#,(s,i,u) + c^ln 





)] 







a: a:' 



{s,t,u,/i h ). (76) 



Note that the hard function is Hermetian, by definition, but not real. Its Hermeticity is 
preserved by the renormalization group evolution. 

In the remainder of this section, we give the explicit expressions for ch, Jh, Mu, and A^ 
for the various channels. 



5.1 qq — >• qq channels 

The values of ch and 7^ for all of the qq — > qq crossings are 



7* = ) (-12CW + ■ • • , ( 7 



c H = 4C F (77) 

which agree with Eqs. (1681) and (ITUl) . The mixing matrix for the representative channel qq' — > 
qq' is 

/4C F [£(u)-£(»)]-C A [L(t) + L(«)-2L(s)] 2[L(«) - L(s)]\ 

M " (s ' i,ll) = l 8fW«)-L(.)] J- (79) 

and the mixing matrix for other channels can be obtained by using the crossing relations given 
in Table [TJ 

For conciseness, since the matrix does not depend on /i, we will use the abbreviated nota- 
tion, following [19], 

T = L(t) - L(s) (80) 



19 



U = L(u)-L{s) . 

Then 

Mu(8,t,u) = I ^ 

When identical particles are involved, one has to be careful to use the RGE appropriate 
for that channel. Recall that we used a Fierz identity to rewrite the O sut operators as linear 
combinations of the O stu operators in Eq. (I19p . The operators O sut and O stu evolve differently 
and to ensure this linear combination is not altered by evolution, we must have 

M(s, t, u) = B • M(s, u, t) ■ B , (83) 

which is easily seen to hold. In fact, we have already used this relation to write all of the 
contributions to the qq — > qq channel in terms of a basis which evolves with Mjj(s, u,t), thus 
letting us calculate a single hard function for the qq — > qq channel. 
The eigenvalues of Mjj are 



A± = %t/ - T) - ^f±^jUT + \C\{T - Uf . (84) 

Here, ± is shorthand for the two complex roots of the discriminant. The diagonal basis is 
related to the original basis by 

Cl = \ ± C{ + ^UCl. (85) 

The normalizations are arbitrary. To change to the diagonal basis for the Wilson coefficients 
or the hard functions, it is convenient to use matrix notation 

C r ± = F ±I Cj r = (F-C r ) ± (86) 
H ±± , = F ±I Fl,j H U = (F-H- Ft) ±± , , (87) 

with 

F±i(s,t,u) = 

The same matrix applies to all channels, with the crossing given in Table [TJ Note that the 
explicit factor of In ^ = L(t) in Eq. fl67|) must be crossed as well. If t crosses into s, this 
logarithm has an imaginary part given by L(s) in Eq. (I38p . However, this particular imaginary 
part drops out when C l + and are combined into the hard function so it can generally be 
ignored for physical applications. The imaginary parts in F±j and X± are important and must 
be included. 
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5.2 gg — > qq channels 

For the gg — >• qq channels, 



c H = 2C F + 2C A 



and 

Mu(s,t,u) 
The eigenvalues are 
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where Xk are solutions of 

X 3 K - C A {T + U)X 2 K + [C\TU — (T — f/) 2 ] A x + C A (T - U) 2 {T + U) = 



(92) 



(93) 



It is straightforward to solve this this cubic equation, with general complex coefficients, but 
the explicit solutions are unilluminating. The Wilson coefficients and hard function in the 
diagonal basis are 
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(95) 



5.3 gg — > gg channel 

For gg -)■ gg, 



ch 



±C A 

7H = (^) (-4A) + 



(96) 
(97) 



The mixing matrix is given by Kidonakis et al. in [19]. Our color basis agrees with theirs, but 
they have tow compared to the convention of [28] which we also use. In addition, we use 
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log -| rather than log -% for the /i-dependence piece, which changes the diagonal elements of 
the matrix. In our notation, the mixing matrix is 

Mu{s,t,u) = 

fc A {U-2T) 
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-2C A T 




The eigenvalues and eigenvectors of this matrix are given for general N in [19]. For N 
the results simplify, and we find 



A 



The matrix of eigenvectors is 



\ 1 = A 2 = -3T 


(99) 


A 3 = A 4 = 3(U - T) 


(100) 


A 5 = A 6 = 3(U - 2T) 


(101) 


A 7 = 2(U - 2T) 


(102) 


8j9 = 2U - 4T±WT 2 -TU + U 2 . 


(103) 


F KI (s,t,u) = 


(104) 



( 2U 


7T 


2J7 





1 





2U 


-2T 


2U 











7U 


2T 


7C/ 


1 





-1 


1 


1 


1 



2(T - C/) 



2(T - C7) 



7T 

0-10 
7(U-T) -2T 7(U-T) 
1 



ai(A 8 ) a 2 (A 8 ) a 3 (X 8 ] 
^ai(Ag) a 2 (A 9 ) a 3 (A 9 ; 



2(Z7-T) 

1 

a 4 (A 8 ) 
a 4 (A 9 ) 



2T 

1 

^(As! 
05(A 9 ; 



2(rj 






-QU 
1 
- T) 6(T - [/) 



-6U Q(U-T) \ 




1 

a 6 (A 8 ) 
ae(A 9 ) 





-1 

a 7 (A 8 ) 
a 7 (A 9 ) 







-6T 


-1 

a 8 (A 8 ) 
a 8 (A 9 ) 





6T 




-1 

ag(A 8 ) 
ag(A 9 ) 



(105) 



with 

ai(A) 



AT 2 U 2 {T - Uf {AT 2 - ATU + 5f7 2 ) (5T 2 - 5TU + 8f/ 2 ) 
240T 4 + 30T 3 (A - 1QU) + 9T 2 f/(56f/ - 5A) + TU 2 (35X - 2Q4U) + 10U 3 (QU - A) 
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a 2 (A) = -T 2 (T - U) (T 2 (6A - 8U) + T[/(8£7 - A) + 5U 2 (\ - AU)) (106) 
o 

a 3 (A) = ai(A) 

a 4 (A) = \t(T - Uf (T 2 {8U + 6A) - TU{W + 11A) + \W 2 (2U + A)) 

a 5 (A) = a 2 (A) 
ag(A) = 04(A) 

a 7 (A) = (T — £/) 2 (-10T 3 A + T 2 U(1QU + 23A) - TU 2 {2W + 23A) + 10£/ 3 (2£7 + A)) 
a 8 (A) = T 2 (-10T 3 A + T 2 U(1QU + 7A) - TU 2 (6U + 7A) + 10?7 4 ) 
a 9 (A) = 2T 2 (T - J/) 2 (AT 2 - ATU + 5U 2 ) . 

6 Soft function 

The results for the NLO Wilson coefficients and their RG evolution apply for any physical 
process for which massless 2 — 2 scattering gives the leading order contribution. If we are 
interested in a particular dijet observable, we must in addition calculate a process-dependent 
soft function. In general, soft function calculations may be much more difficult than the hard 
function calculation, since they can involve arbitrarily complicated phase space cuts. However, 
as we will show in this section, there are some universal features which must hold for any soft 
function. In particular, the //-dependence of the soft function and its RGE to 2-loop order are 
almost completely fixed. 

In SCET, a soft function is calculated by taking matrix elements of the products of Wilson 
lines which appear in the SCET operators. For example, for the qq' — > qq' operators, these 
are 

Wi = T{(y 4 t r a y 2 )(y 3 W 1 )} 

w 2 = T{(y 4 t iy 2 )(y 3 t iy 1 )} . (107) 

The soft function is then calculated from 

Sjj ({*}, nf ,fi) = J2 (°l ^ W (X,\ m |0> F s ({k}) , (108) 

where the sum is over soft radiation in the final state. The function Fs({k}) encodes the 
dependence on various projections on the soft momenta related to the definition of the ob- 
servable. No matter what the observable is, the soft function can only depend on directions 
nf of the various Wilson lines, and on arbitrary soft scales {k} relevant to the projections. 
Because of factorization, it cannot depend on the energy of the jets, the hard scales s, t, u or 
the energy fractions Xi of the PDFs. 
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For the final cross section to be independent of /i, the soft function must satisfy a renor- 
malization group equation closely related to that of the hard Wilson coefficients. The soft 
function evolution will not be local in {k}, but must be local in Laplace space for factorization 
to hold. The Laplace transform of the soft function is defined, for one scale k, as 

Su(Q,n^,fx) = J d&exp (-JL-j Sufcn?,!!) , (109) 

g with the natural generalization to more scales {k}. In Laplace space, the jet function and 
PDF evolution is also local. Furthermore, since the jet functions and the PDFs evolve diago- 
nally in color, the color mixing terms in the soft function evolution must exactly compensate 
the color mixing terms in the hard function evolution. 
It follows that the soft function must satisfy 



d 



dln/i 

where 



S u ({Q}, nl /i) = -S IL ({Q}, <, //) T S LJ - FflSlj ({Q}, nj\ /*) , (110) 



rfj = ^7cus P Cq In + 7 cusp r(nf ) + 75^ S u + 7cuspM 7 j(nf ) . (Ill) 

This has many of the same properties as the hard function RGE, Eq. fl67|) : 

1. The //-dependence is proportional to 7 CUS p times the identity matrix in color space. 

2. The color mixing is proportional to 7 CUS p and fixed by a single matrix. 

3. The unknown quantity 75 depends only on the channel, not the color structure. 
Furthermore, for the soft function we can observe further that 

4. The matrix Mu{nj) can be taken to be the same as the Mjj matrix appearing in 
Eq. (157)1 . which depends on s,t and u. This is possible because only dimensionless 
Lorentz-invariant ratios appear, and (rij • rij)/(nk -ni) = (pi • Pj) / (pk • Pi) ■ To make this 
transparent, we can even write 

M = " E ^Y 1 [ L ^ - L M (112) 

where = ±rij -rij, taking the + sign when rii and rij are both incoming or both 
outgoing, and the — sign otherwise. This form of M is a variation of the form presented 
in [201. 



5. The color-diagonal piece, proportional to 7 CUS p^/j, m ay have dependence on the directions 

and on the soft scales {Q}- 

6. RG invariance lets us solve for 75 in terms of 7^ and the jet and PDF anomalous 
dimensions 7j and 7/. 
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Note that in denning the matrix Mjj with the hard function, Eq. ( 16 7p . there was an 
implicit convention in the diagonal terms. For example, if we had written In ^ instead of 
In =| the diagonal terms in Mjj would be different. This change of convention can then be 
absorbed in the r(n^) term in the soft function anomalous dimension. The terms in the 
soft function evolution equation proportional to 5jj are also different from those in the hard 
function evolution equation because they must additionally compensate for the evolution of the 
jet functions and PDFs. Since the number and type of jets which contribute depends on the 
observable, this term is observable-dependent. Once an observable is defined, the entire soft 
function RGE can be fixed either by a 1-loop computation or by analyzing the factorization 
formula. 

Since the mixing terms in the soft function are proportional to Mjj, there is no mixing in 
the basis in which Mjj is diagonal. In the previous section, explicit expressions were given for 
the matrices Fjk which let us change from the / basis, in which the NLO matching is simple, 
to the K basis, in which the evolution is simple. We had written 

H KK , = (F ■ H ■ F^) KK , . (113) 

Consequently the soft function transforms as 

S KKI = [{F-y -S-F- l ] KK ,. (114) 

The cross section involves a sum over color structures 

Y^HijSji = J2 Hkk ' Sk ' k = T±( H S) , (H5) 

IJ KK' 

which is basis independent. Since both H and S are Hermetian, this sum is real. However, 
note that F may not be unitary since Mjj is in general not Hermetian. 
The general solution to the soft function RGE in the diagonal basis is 



S KK/ ({Q},nf,iS) = exp - 2c Q S(fi s ,fx) + 2A s (fJ. a , fj) 



x exp 



A r (fis, fJ.)[\ K + \* K > + r(nf) + r«)* + 2c Q In 



{0} 

Us 



Skk>({Q}X,Hs)- (116) 



This can be transformed back to momentum space using the techniques described, for example, 
in [10], but the general expression is not particularly illuminating. Note that even in the K 
basis, in which the evolution is diagonal, the soft function itself may still have non-diagonal 
terms (K ^ K') thus leading to color mixing among the Wilson coefficients C V K . For next-to- 
leading-log resummation (NLL), the soft function is completely determined by its anomalous 
dimensions, and then Skk'{{Q}, /-0 is diagonal for any /x. However, at NNLL the additional 
non-diagonal terms in the finite parts of the soft function have to be included. Nevertheless, 
due to the very general arguments presented above, each component of the soft function at the 
scale \i is proportional to the the same component at any other scale //, Skk'(^) oc Skk'^h'), 
as this equation shows. 
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Although a factorization formula will let us deduce 75 from RG invariance, it is satisfying 
to also calculate this anomalous dimension directly. Such a calculation provides an important 
check that the observable is well defined and that the factorization theorem actually works. 
Moreover, performing the calculation of a particular soft function to order a s will also produce 
the //-independent a s terms, which are necessary for NNLL resummation. An example can be 
found in [38]. Nevertheless, it is easy to imagine a complicated enough observable for which 
direct calculation is not possible. In that case, these general results will allow us to resum the 
next-to-leading logarithms (NLL) without actually doing the soft integrals. 

7 Conclusions 

We have presented the tools necessary for an effective field theory calculation at NNLL of 
observables based on 2 — > 2 processes in QCD. The 1-loop matching coefficients for all the 
2—7-2 processes have been compiled in a format that is conducive to the computation of 
dijet observables in SCET. The evolution of these matching coefficients is independent of the 
observable and we have given analytic formulas for the Wilson coefficients at an arbitrary scale 
/1 including NNLL resummation. 

For each channel, there is a basis in which the evolution of the Wilson coefficients is 
diagonal in color space. These bases and the associated eigenvalues of the evolution equations 
were given explicitly for all channels. Unfortunately, the basis in which the 1-loop matching is 
simple and the basis in which the renormalization group evolution is simple do not coincide. 
Nevertheless, since we have provided the explicit matrices to go between bases, numerical 
evaluation of the resummed expressions is now straightforward. 

General soft functions, which may depend on the observable, were also examined. They 
were shown to have similar universal properties to the hard function as required by the consis- 
tency of the effective field theory. The bases which diagonalize the Wilson coefficient evolution 
also diagonalize the soft functions evolution and simplify the form of the resummed expres- 
sions. 

With these results, it will be possible to begin accurate calculations of dijet observables 
including NNLL resummation. As a start, it would be useful to compute a simple observable 
which can check the RG evolution for the soft function we have derived here through explicit 
calculation. Such a calculation is given in [38]. Then, more interesting observables may be 
approached, including, for example, the jet px spectrum at large pt, jet substructure, such as 
jet masses or angularities [31], hadronic event shapes [39], or even possibly observables related 
to the color flow in an event [3D] . 
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